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Following the early discovery of stripe-like order in La-based copper-oxide 
superconductors, charge ordering instabilities were observed in all cuprate 
families. However, it has proven difficult to distinguish between uni- (stripes) 
and bi-directional (checkerboard) charge order in Y- and Bi-based materials. 

Here we use resonant x-ray scattering (RXS) to measure the two-dimensional 
structure factor in YBa2Cu306+y, in reciprocal space. Our data reveal the 
presence of charge stripe order, i.e. locally unidirectional density waves, sug¬ 
gesting it as the true microscopic nature of charge modulations in cuprates. At 
the same time, we find that the well-established competition between charge 
order and superconductivity is stronger for charge correlations across than 
along the stripes, which provides additional evidence for the intrinsic unidi¬ 
rectional nature of the charge order. 

Recent studies of Y-based copper oxides have highlighted the importance of a charge- 
ordered electronic ground state, also termed charge-density-wave (CDW), as a central ele¬ 
ment within the phenomenology of high-temperature superconductors { 1 - 14 ). The family of 
YBa2Cu306+y (YBCO) compounds have yielded a wealth of experimental results that enabled 
advancements in our understanding of CDW instabilities and their interplay with superconduc¬ 
tivity { 9 - 11 , 15 - 22 ). 

YBCO is a layered copper-oxide-based material where hole doping is controlled by the oxy¬ 
gen stoichiometry in the chain layer - characterized by uniaxial CuO chains running along the 
crystallographic b axis. In addition to ordering within the chain layer - attained via the periodic 
alternation of fully-oxygenated and fully-depleted CuO chains - recent experiments have ex¬ 
tensively shown the presence of charge ordering in the Cu 02 planes, with an incommensurate 
wavevector (5~0.31 reciprocal lattice units { 23 ), corresponding to a period of approximately 3 
unit cells in real space { 9 - 11 ). Although the stripy nature of La-based cuprates has been long 
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established (1-3), the loeal symmetry of the CDW in YBCO has not yet been resolved as both 
eharge stripes (in the presenee of 90° rotated domains)and eheekerboard are eonsistent with the 
globally bidireetional strueture of the CDW - eharaeterized by waveveetors along both the a 
and b axes, at Qa ~ (0.31,0) and Qh ~ (0,0.31), respeetively (10,11,24-26). This leaves open 
the fundamental question as to whether stripes are the underlying eharge instability in the whole 
elass of hole-doped euprates. 

Here we study the loeal density eorrelations of the eharge-ordered state, and the interaetion 
of the latter with supereonduetivity (SC) in underdoped YBCO, using resonant x-ray seatter- 
ing (RXS). This teehnique, whieh is now at full maturity, represents a unique eombination of 
diffraetion (to probe reeiproeal spaee) and resonant absorption (allowing element-speeiheity 
and therefore site-selectivity), and directly measures the structure factor S(Qx, Qy), where 
and Qy represent the momenta along the reciprocal axes H and K, respectively. The structure 
factor is linked to the density-density correlation function, and therefore to the CDW order 
parameter in momentum space (see Materials and Methods in Ref. 27). Here, in order to recon¬ 
struct the two-dimensional (2D) structure factor with high resolution, we utilize a specifically 
devised RXS probing scheme whereby a charge-ordering peak is sliced along different direc¬ 
tions, parametrized by the azimuthal angle a (Fig.[T]). The resulting 2D shape of the CDW peaks 
rules out checkerboard order and is instead consistent with a stripy nature of charge modulations 
in YBCO (28). We carry out RXS measurements around the CDW waveveetors ~ (0.31, 0) 
and Q. (0,0.31) for three detwinned, oxygen-ordered YBCO samples - YBa2Cu306.5i 
(Y651, with hole doping p ~ 0.10), YBaaCugOg.er (Y667, p ~ 0.12), YBaaCugOg.Ts (Y675, 
p ~ 0.13). In our experimental scheme the CDW peaks are scanned in a radial geometry via 
control of the azimuthal angle a (29) (Fig.[^). At the Cu-L^ edge, the measured signal is 
mainly sensitive to periodic variations in the Cu-2p —)■ 3d transition energy (30, 31), i.e. to 
a scalar quantity, even though the detailed contribution of a pure charge modulation vs. ionic 
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displacements to the RXS signal eannot be deeoupled (32). In addition, the poor eoherenee of 
the CDW aeross the Cu02 planes (JJ, 24, 26) qualifies this eleetronie ordering as an exquisite 
two-dimensional phenomenon, thus motivating our foeus on the strueture faetor in the (Qx^ Qy) 
plane. Representative seans of the CDW peak for different a values and at the supereondueting 
eritieal temperature T ~ Tj. are shown in the insets to Fig.[^A for the Q;, and Qq CDW peaks 
of a Y667 sample, respeetively. A ehange in the peak half-width-at-half-maximum (HWHM) 
AQ between a = 0° and 90° is already apparent, but is even better visualized in the eolor map 
of Fig. whieh shows the sequenee of Q-seans vs. azimuthal angle and the eorresponding 
variation of AQ for in the range a = —90° to 90°. This same proeedure is applied to all 
three YBCO doping levels, for both the Qa and Qb CDW peaks; polar plots of AQ vs. a are 
shown in Figs.[^C-E for Y651, Y667 and Y675, respeetively. With the aid of the ellipse fits 
to the CDW profiles (eontinuous lines), four key aspeets of these data stand out: (i) all peaks 
show a elear anisotropy between the two perpendieular direetions a = 0° and 90°; (ii) for eaeh 
doping, the and Q;, peaks have different shapes, and in the ease of Y651 and Y667 this is 
even more evident as the peaks are elongated along two different direetions; (hi) the departure 
from an isotropie ease, quantified by the elongation of the CDW ellipsoids, inereases towards 
the underdoped regime and is opposite to the evolution of orthorhombieity, whieh is instead 
maximized at optimal doping (see Supplementary Materials and Fig. S4 for a more detailed 
diseussion); (iv) the peak elongation at and Qb, evolving from biaxial (Y651 and Y667) to 
uniaxial (Y675), is ineonsistent with the doping-independenee of the uniaxial symmetry of the 
CuO ehain layer, whieh rules out the possibility that the observed CDW peak strueture is ex- 
elusively eontrolled by the erystal’s orthorhombie strueture - however, the uniaxial anisotropy 
observed for Y675 might refleet a more pronouneed interaetion between the Cu-0 planes and 
ehains, possibly due to the inerease in orthorhombieity upon hole doping (see more detailed 
diseussion in the Supplementary Materials). 
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The observed 2D peak strueture direetly bears the signatures of a breaking of fourfold (C4) 
symmetry at both the macro- and nano-scale, and therefore of the emergenee of a stripe-ordered 
state. In faet, under the fourfold symmetry the eleetronie density would be invariant under a 90° 
rotation in real spaee (x —)■ y, y —)■ —x), whieh is equivalent to a 90° rotation in momentum 
spaee {Q^ —)■ Qy, Qy — )■ —Qx)- Instead, the CDW strueture faetor S{Qx,Qy) is elearly not 
invariant under sueh operation, as shown in the bottom-right diagrams in Figs.[2 C-E whieh 
eompare the original S{Qx, Qy) to their 90° rotated version S{Qy, —Qx)- This finding demon¬ 
strates an unambiguous breaking of global symmetry in all investigated samples, and might 
elueidate the origin of the anisotropy observed in the Nernst effeet (20) and in optieal birifren- 
genee measurements {33). The real-spaee representation of eharge order branehes off into two 
possible seenarios: (i) a biaxial anisotropy, where both x- and y-elongated domains {34) are 
present (Figs.[^ A and C); (ii) a uniaxial anisotropy, where only y-elongated (or, equivalently, 
x-elongated) domains are found (Figs.|^ B and D). Note that these domains need not neeessar- 
ily lie in the very same layer, but they need to be present at the same time within the bulk of the 
material (e.g., they ean be present in alternating layers, while still leading to the same momen¬ 
tum strueture). The momentum-spaee representation of the order parameter - and therefore of 
the eleetronie density fluetuations - is shown in the eorresponding panels in Figs.|^ E-H. Here 
S'a(Q) (red ellipses) and S'6(Q) (blue ellipses) represent the strueture faetor assoeiated to eharge 
modulations along a and b, respeetively. The profile of a single strueture-faetor peak is eom- 
pounded with two eontributions: the underlying CDW symmetry as well as its 2D eorrelation 
length, whieh ean also be anisotropie. As a net result, the anisotropy of a single peak in Q-spaee 
eannot be used to diseriminate between different CDW symmetries. Instead the latter ean be 
resolved by probing the 2D CDW strueture faetor, i.e. by eomparing the CDW peak shape for 
both Qa and Q;,. By inspeeting the diagrams in Figs.[^ E-H, one ean reeognize a eommon 
trait of eheekerboard struetures in momentum spaee, in that the following eonditions (Figs.[^ 
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G and H) must always be verified by symmetry: A(5“ = and AQ“ = AQy - i.e. the peak 
broadening along a given direction must coincide for A^ and Ab (see bottom of Figs.E-H for 
case-specific conditions on the peak linewidths). Intuitively, this follows from the fact that - for 
the checkerboard case - the charge modulations along a and b axes must be subject to the same 
correlation lengths within each domain - irrespective of its orientation - and therefore lead to 
an equivalent peak broadening along the same direction in reciprocal space, in contrast to our 
findings for the CDW linewidths (see Supplementary Materials for more details). From this 
symmetry analysis, we can conclude that for both uniaxial and biaxial anisotropy it is in prin¬ 
ciple possible (35) to discriminate between a pure checkerboard and a pure stripe charge order, 
even in presence of a distribution of canted domains (see Tables SI-S3 for a complete classifi¬ 
cation). Ultimately, the inequivalence of the peak broadening AQ along different directions for 
all studied YBCO samples, combined with the macroscopic C 4 symmetry breaking, provide a 
clear cut evidence for the unidirectional (stripe) intrinsic nature of the charge order (28). 

Having established the underlying stripe-like character of charge modulations in the CUO 2 
planes, we turn to the temperature dependence of the longitudinal and transverse correlation 
lengths, respectively parallel and perpendicular to the specific ordering wavevector. These can 
be extrapolated by inverting the momentum HWHM AQ, as illustrated in Fig.|^. Fongitudinal 
(transverse) correlations are then given by = AQu"^ (^_l = AQj^), where AQu (AQ^) repre¬ 
sents the momentum linewidth in the direction parallel (perpendicular) to the ordering wavevec¬ 
tor. We subsequently study the temperature dependence of .^n and .^_l for both the Qa and Qb 
ordering wavevectors (Figs.[^ B-F). First we observe the rise of correlation lengths below the 
CDW onset near 150 K, followed by their suppression below the SC transition temperature Tc, 
which confirms the competition between these two orders, in agreement with recent energy- 
integrated as well as energy-resolved RXS studies (10,11,24,25,30). However, the drop in the 
correlation lengths below Tc (A0 is in all instances larger for the longidutinal correlations, or 
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A?ll > In particular, the discrepancy between and A^^, while small for Y675, is 

quite substantial for the more underdoped Y667 and Y651. This anisotropy provides additional 
evidence for the unidirectional nature of the charge ordering and thus the breaking of C 4 sym¬ 
metry, since a bidirectional order would exhibit an isotropic drop in correlation length across T^. 
This anisotropy has an opposite doping trend from the crystal orthorhombicity whose associated 
anomalies across Tc increase with hole doping [as revealed e.g. by lattice espansivity measure¬ 
ments {36)]. The inferred real-space representation of the evolution across Tc is schematically 
illustrated in Figs.|^ G and H, where nanodomains are used to pictorially represent a charge- 
ordered state with finite correlation lengths. We conclude that the largest change occurs along 
the direction perpendicular to the stripes. This reflects the tendency of the SC order to gain 
strength, as temperature is lowered, primarily at the expense of longitudinal CDW correlations, 
suggesting that the mechanism responsible for the density fluctuations across the periodically 
modulated stripes might be the main one competing with the Cooper pairing process. 

Our results may explain many common aspects between CDW physics in YBCO and the 
stripy cuprates from the La-based family, such as: thermoelectric transport {37, 38)-, strength 
of the order parameter {39)\ out-of-equilibrium response {40, 41); energy-dependent RXS re¬ 
sponse {30, 31). The nanoscopic nature of the stripe instability and the presence of both a- and 
b-oriented domains also clarify why this broken symmetry has been difficult to disentangle from 
a native bi-directional order {10,11,42), therefore requiring a tailored experimental scheme to 
resolve the 2D CDW structure factor S{Qx,Qy)- The pronounced directionality in the com¬ 
petition between superconductivity and stripe order reveals new key aspects of the underlying 
interplay between particle-particle and particle-hole pairing in high-temperature superconduc¬ 
tors, and provides novel insights for an ultimate understanding of these materials. 
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Figure 1: Charge order topology in momentum space. (A) Schematic representation of 
the momentum structure of charge modulations in YBCO. Left inset: Selected momentum 
scans of the CDW peak along the b axis at Q;, = (0,0.31), for different azimuthal an¬ 
gles (a = 0°,45°,90°); continuous lines represent Lorentzian fits, horizontal bars denote the 
linewidth /S.Q (HWHM). Right inset: Same as for the left inset, but for the CDW peak along 
the a axis at = (0.31,0). (B) Color map of a series of Q-scans (normalized to the peak 
height) slicing the Qb peak between a = —90° and 90°; black bars represent the linewidth 
AQ, which is largest at a = 0°. (C-E) Polar plots of /\Q as a function of the azimuthal angle a 
for Qa (red) and Qb (blue) in YBa2Cu306.5i (Y651), YBa2Cu3 06.67 (Y667), and YBa2Cu306.75 
(Y675), respectively. Concentric grey circles are spaced by O.OlA ; continuous lines are fits to 
an elliptic profile. Bottom-right insets: CDW peaks represented as filled ellipses and compared 
with their rotated versions (hollow ellipses), for each doping. 
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Figure 2: Domain mesostructure in real and reciprocal space. (A,B) Stripy domains along a 
(red stripes) and b (blue stripes), in the presence of biaxial (A) and uniaxial (B) anisotropic cor¬ 
relations. The domains are visualized in the same layer, although a situation in which they are 
present in alternating layers in a 90°-rotated arrangement is equally possible. (C,D) Checker¬ 
board domains in the presence of biaxial (C) and uniaxial (D) anisotropic correlations. (E-H) 
Corresponding structure factors in reciprocal (Q) space. In case of the simultaneous presence of 
both CDW components (checkerboard), the imprinted correlations must be equal for the density 
wave along a and b, thus imposing an equivalent peak structure at and (panels G and H). 
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Figure 3: Interplay between stripe order and the superconducting state. (A) Schematics of 
the CDW structure factor in Q-space, illustrating how longitudinal (^y) and transverse (^_l) cor¬ 
relation lengths are extracted from the data by taking the inverse of the peak HWHM AQ. (B-F) 
Temperature dependence of (green markers) and ^_l (orange markers) for the various sam¬ 
ples and CDW peaks investigated (see legend). Continuous lines are guides-to-the-eye; double 
arrows highlight the drop in correlation lengths A^y and below Tc. (G,H) Cartoons rep¬ 
resenting the evolution of the stripy nanodomains from above (G) to below (H) illustrating 
how longitudinal correlations undergo a larger suppression in the presence of superconducting 
order. 
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Materials and Methods 

Sample characterization. This investigation was performed on three detwinned, oxygen- 
ordered underdoped YBa2Cu306+j; single crystals: (i) y = 0.51, p ~ 0.10, Tc = 57 K, Ortho II 
oxygen-ordered, hereafter labeled Y651; (ii) y = 0.67, 0.12, Tc = 65 K, Ortho VIII oxygen- 

ordered, hereafter labeled Y667; (iii) y = 0.75, p ~ 0.13, Tc = 70 K, Ortho III oxygen-ordered, 
hereafter labeled Y675. The in-plane lattice constants, which in turn determine the degree of or- 
thorhombicity, are: (i) a~3.8362A and 6~3.8740A in Y651; (ii) a~3.8257A and 6~3.8814A 
in Y667; and (iii) a ~ 3.825A and b ~ 3.885A in Y675. The superconducting critical tem¬ 
perature Tc was determined from magnetic susceptibility measurements. The in-plane lattice 
constants have been measured using XRD in Y651 and Y667, and are taken from Ref. 44 for 
Y675. Details on the sample growth and the Tc-to-doping correspondence can be found in 
Ref. 45. 

Resonant soft X-ray scattering. Resonant x-ray scattering (RXS) is a coherent probe of fluc¬ 
tuations in those physical quantities (X) that are directly or indirectly coupled to light, and is 
therefore proportional to the Fourier transform of the correlation function of 5X, i.e.: 

Irxs{Q) oc j dre-*Q-' j dr' {5X{r')5X{r' + r)) = (<5X(-Q)<5X(Q)), (SI) 

where Cx(r) = / dr' ((iX(r')(5X(r' -f r)) represents the correlation function of 5X. 

Charge ordering is manifested through variations in different physical quantities, such as the 
ionic displacements (X = (5R, in this case X is a vector operator), the local valence (X = p), 
or the energy shifts of the resonant transition (X = AE). The sensitivity of the experimental 
signal to these fluctuating quantities is material-dependent: in the cuprates it has been shown 
that the resonant scattering is predominantly controlled by variations in the energy shifts (i.e., 
a scalar field) [30,31], whereas in a different material, e.g. IrTe 2 , the RXS signal (at the Te-M 
absorption edge) was shown to be primarily associated to a combination of valence modulations 
and atomic displacements [46]. In presence of valence modulations or energy shifts (i.e., scalar 
operators), the fluctuating quantity is proportional to modulations in the electronic density p, 
i.e. 6X oc 6p, which leads to: 

iRxsiQ) « {Spi-Q)Spm = {6p*iQ)6pm = SiQ), (S2) 

where the second equality follows from the fact that the density operator is a real function, and 
S'(Q) is termed the static structure factor. 

For this work, RXS measurements were performed on a 4-circle diffractometer in a 10“^° mbar 
ultra-high-vacuum chamber, with a photon flux around 5 ■ 10^^ photons/s and ^ ~2 ■ 10“^ en¬ 
ergy resolution. In addition, fully polarized incoming light is used, with two available configu¬ 
rations: a (polarization vector perpendicular to the scattering plane) or tt (polarization vector in 
the scattering plane). In order to maximize the charge order signal, the series of measurements 
as a function of the azimuthal angle a were performed at the peak energy of the Cu-Ls edge 
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(hu = 931.5 eV, see Fig. S[^, at a deteetor angle 9det = 170°, and at a temperature near the su- 
pereondueting transition Tc. All RXS seans were measured with a multi-ehannel-plate deteetor 
with an angular resolution A9 ~ 0.2°, eorresponding to an equivalent momentum resolution 
AQ ~ 0.0023A"^ at hz/ = 931.5 eV. 

One should note that the speetrometer integrates the energy of seattered photons; however, 
sinee inelastie proeess are incoherent in nature, the only portion of the speetrum whieh earries 
the information on the momentum strueture of the eleetronie density is the elastie (zero-loss) 
line. Sinee the inelastie eontributions, albeit strong at or near resonanee, do not have a sharp 
strueture in momentum spaee, the momentum strueture of the CDW peaks as observed with an 
energy-integrated speetrometer is still predominantly determined by the elastie proeesses, and 
therefore it is representative of the statie eleetronie modulations in Fourier spaee. This is also 
eonfirmed by the elose agreement between the CDW momentum linewidth previously measured 
in the energy-resolved (RIXS) and energy-integrated (RXS) mode (see Fig. 4 in Ref. 10). 

Supplementary Text 

Azimuthal scans and experimental geometry. The azimuthal angle a is defined as the angle 
between the direetion of the RXS sean in momentum spaee and the erystallographie a axis, 
or equivalently its reeiproeal axis H. This eonfiguration is illustrated in Fig. S|^, whieh also 
highlights the euts in momentum spaee (for different a values), at the eharge-density-wave 
(CDW) ordering waveveetors studied here, namely = (0.31,0,1.5) and Q;, = (0,0.31,1.5). 
The value of L = 1.5 is ehosen to maximize the CDW intensity [10]. However, as studied in 
detail in reeent work [11,24,26], the dependenee of the CDW ordering peak on the out-of-plane 
waveveetor eomponent L is very weak, a signature of the two-dimensional nature of the eharge 
order. For this reason, only the planar momentum strueture of the CDW order parameter is 
studied here, and we will heneeforth refer exelusively to the in-plane projeetion of the CDW 
peaks, and use the notation Qa = (0.31, 0) and Qfe = (0, 0.31). 

In order to extraet a proper profiling of the CDW peak, it is important that eaeh Q-sean 
outs through the CDW peak maximum. To ensure that this eondition is met for eaeh azimuthal 
angle a, angular optimization perpendioular to the Q-sean direetion was performed by soanning 
through the peak as a funotion of the transverse angle x (see lower-left inset of Fig. S|^). By 
finding the x angle maximizing the CDW signal, we have been able to eut preeisely through the 
peak maximum for every azimuthal position. 

In a RXS measurement, the angle-to-momentum oonversion is enooded in the equations: 



(S3) 


where Q is the waveveetor magnitude, Q\\ is its projeetion on the sample surfaee (a — b plane), 
while 6*det is the deteetor angle, 6'sampie is the sample angular position (see sehematios of the 
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probing geometry in Fig.|^), while u is the photon frequeney. For the eharge ordering peak, 
the angular values are the following: 6*det = 170°, 6*gampie = 32.5°. 

The resolution of the RXS diffraetometer is evaluated by measuring the Bragg refleetion 
(0,1,3) at high photon energy {hu = 2230 eV). By seanning the deteetor arm aeross the Bragg 
refleetion - whose intrinsie broadening is mueh smaller than the instrumental resolution and 
therefore represents a 5-like signal - we estimate the angular resolution to be A6*det ~ 0.2° (see 
lower part of Fig. S|^). In order to estimate the momentum resolution, we have eonverted the 
angular resolution sean (see again bottom part of Fig. S|^) into the surfaee projeeted wavevee- 
tor and then evaluated AQ|| using Eqs.lSSjas follows: 


AQll Q\\ (^det T A0det/2, ^sample) hu'j Q|| (^det A0det/2, ^sample; 

= Q|| (170° + 0.1°, 32.5°, 930.5el/) - Qy (170° - 0.1°, 32.5°, 930.5el/) 

Att 

= — [sin (85.05°) cos (180 - 170.1 + 32.5°) - sin (84.95°) cos (180 - 169.9 + 32.5°)] 

Ch> 

= 0.694668 - 0.692344 ~ 0.0023A“\ (S4) 


The value obtained amounts to approximately the 5% of the average full-width-at-half-maximum 
for the eharge order peak (whieh is of the order of 0.04 — 0.05A ). Sinee the two eontributions 
add in quadrature, if we use the term to denote the intrinsie eontribution as opposed 

to the one from resolution, then we have for the measured linewidths that AQ^^p = 

^QLtive + Since AQ^^, ~ 0.05 ■ AQ^,,tive> we obtain that = 0.9987 ■ AQ^^p, 

which shows that approximately the 99.8% of the observed peak broadening (see again Fig.l of 
the main text) represents intrinsic sample properties. 


Data analysis methodology. RXS scans have been analyzed by means of a least-square non¬ 
linear regression analysis, with fitting function defined as the sum of a Lorentzian peak and a 
cubic background polynomial function B\ 


7rxs — 


A 


Q-Q^ 

AQ 


+ B{Q) 


+ 1 


B{Q) — do + Q + Q^i 


(S5) 


where Q* is the peak position, A is the amplitude and AQ is the half-width-at-half-maximum. 

All momentum scans were measured around the ordering wavevector using two light polar¬ 
izations, vertical (a) and horizontal (vr). A preliminary analysis has shown that the linewidths do 
not depend on the light polarization, consistently with the fact that the linewidth is an intrinsic 
property of the system (related to correlation lengths, as will be discussed later) and therefore 
does not vary upon changing an external parameter, such as light polarization. In order to pro¬ 
vide a more precise quantitative estimate of the linewidth AQ with a reduced statistical error, 
this parameter has been subsequently constrained to be the same during least-squares fitting of 
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the a- and vr-polarized scans. An example of the fitting analysis is provided in Fig. S|^ 


Stripe order with anisotropic correlations. As anticipated in the Materials and Methods 
section, in our case the RXS observable is proportional to the static structure factor S'(Q): 

/rxs(Q)oc^(Q) = (|<5p(Q)0. (S 6 ) 

Here we focus on how a rotational symmetry breaking in real space - a tell-tale of a stripe¬ 
like charge-ordered state - can be equivalently detected by looking at the structure of density 
modulations in reciprocal space. Whenever (^4 symmetry is preserved, the electronic density 
fluctuations 6 p must be invariant under a 90° rotation in real space (defined through the trans¬ 
formations X ^ y and y —>■ —x), or 6 p{x,y) = 6 p{y, —x). This can be readily translated in 
momentum space by taking the Fourier transform of Sp{x, y), which yields: 


Sp{Qcc,Qy) = j 

= y dx dy -x) 

= j dx' dy' e-^^-^^y'+^y^'^ 6 p{x\y') 
^P{Qyy Qx\ 


(S7) 


which proves that the structure factor S'(Q) has to likewise be invariant under a 90° rotation 
in momentum space, or S{Qx,Qy) = S{Qy, —Qx)- We then apply this procedure to the full 
2D momentum structure of the CDW structure factor as measured using RXS (see Fig. 1 in the 
main text). 

The resulting peak structure is readily shown to violate (^4 symmetry, as illustrated in the 
three panels of Fig. S|^ (reproduced from Fig. 1 of the main text). The extent of the rotational 
symmetry breaking is here quantified via a C 4 , symmetry breaking parameter Ac' 4 , defined as 


Ac4 = 


^(Q,)U^(Qfe)-^(Q,)n^(Qb) / ^(QJU^(Qfe) 


, or equivalently as the area of 


the union of the two CDW peaks (one of which rotated) S{Qa) U S{Qb) rninus their overlap 
area S'(Qa) fl ^(Q;,), divided by their union. The doping dependence of is then shown in 
Fig. S|^. On the same plot, we report the values of the peak anisotropy for both the Qa and 
Qb CDW peaks, as derived by taking the ratio between the peak width along the parallel and 
perpendicular direction with respect to the ordering wavevector, or 7 = AQ||/AQj^. A value of 
7 = 1 yields isotropic peaks, while the larger the deviation from 1 , the more anisotropic the peak 
is (in one or the other direction). Interestingly, for Y651 and Y667, both 7 ^ and % are smaller 
than 1, corresponding to CDW peak structures elongated differently between the two directions 
in reciprocal space. On the other hand, in Y675 we find that 7 a > 1, 7 ?, < 1, corresponding to 
having Qa and Qb elongated along the same direction (see again Fig. S|^). Herefater we will 
use the denomination biaxial anisotropy for the case of Y651 and Y667, and uniaxial anisotropy 
for the case of Y675 (see more below on uni vs. biaxial). 
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More detailed doping-, azimuthal angle-, and temperature-dependent results on peak widths 
and eorrelation lengths are presented in the next two seetions. 


Azimuthal-dependent experimental data and parameters. The two-dimensional strueture 
faetor near the ordering waveveetors Qa and Q;, is well-approximated by a two-dimensional 
Lorentzian funetion defined as follows: 


t2D 

-^Lor 


(Q X') Qy 


A 

( Q.-Ql V , [ Qy-Qi y . . 
\ AQ^ J AQ^ J ^ 



A 



+ 1 


(S8) 


where Q* = (Q*, Q*) is the CDW waveveetor, AQ^ y are the half-widths-at-half-maximum 
(equal to the inverse eorrelation lengths) along x and y, respeetively, while in the last equality 
the redueed eoordinate Q' = Q — Q* was used. By parametrizing Q' in polar eoordinates as 
{Q'x, Q'y) = Q'(cos a, sin a), we ean subsequently study the direetional dependenee of the RXS 


intensity from Eq. S8 where the direetion in reeiproeal spaee is eontrolled by the angle a: 


Ilor (Q, «) = 


A 


A 


Q 


a 


cos^a 




+ 


sm^a 




+ 1 


Q 


,/2 


(S9) 


AQ{ 


+ 1 


Eq.|^is used to fit the azimuthal series of RXS seans, whieh allows extraeting the experimental 
a-dependent linewidths AQ^ (error bars are eomputed from a least-squares fitting analysis), 
shown in Eig. S[^ The same equation also allows reformulating AQ^ as a funetion of the two 
independent parameters AQ^ and AQy, whieh eorrespond to the linewidth at a = 0° and 90° 
degrees, respeetively: 


AQ„ 


cos^a 


AQ 


+ 


sin^a 


-1 




(SIO) 


This last formula is used to fit the a-dependent linewidths displayed in Eig. S|^ the eorrespond- 
ing fit profiles are overlaid as eontinuous lines. The values of AQ^ (AQy) are reported beside 
the horizontal (vertieal) bars in the top-right insets, whieh also define the peak shape (full ellip¬ 
soid) and eorresponding Qx/Qy anisotropy. 

Erom the values of AQ^ we ean derive the a-dependent eorrelation lengths upon simple 
inversion: = AQ~^. The eorresponding data points are plotted in Eig. S[^ Here we also 

introduee longitudinal (.^y) and transverse (.^x) eorrelations, representing the eorrelation lengths 
along the two axes that are respeetively parallel and perpendieular to the speeifie waveveetor. 


Analysis and categorization of possible domain structures. A more detailed study of the 
native anisotropy of the peak shape and eonsequently of the eorrelation lengths, in relation 
to possible CDW domain struetures, is also important in revealing the stripy nature of eharge 
ordering in YBCO. Sueh analysis relies on a few preliminary observations: 
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• The x-ray beam spot size (d ~ SOO/rm) is mueh larger than the average eorrelation length 

~ 40 — 60A), henee variations in the latter with the azimuthal angle a eannot be 
attributed to the (minimal) motion of the beam spot on the probed region of the sample. 

• The shortest and longest eorrelations always oeeur at either a = 0 (eorresponding to 
a direetion parallel to the a axis) or a = 90 (parallel to the b axis) and therefore rule 
out the influenee of extrinsie effeets related to the probing geometry in determining the 
a-dependenee of the eorrelation lengths - whieh would not preferentially seleet a high- 
symmetry direetion, and would rather eause the signal to be minimized or maximized at 
a random angular position. 

• The different variation of ^ with azimuthal angle between Y675 and the other two samples 
(see again Fig. S[^ also exeludes any systematie effeet to the probing geometry. 

Based on the results shown in Fig. S[^ anisotropie eorrelations are eonfirmed to be present 
in all YBCO samples examined. In general, we note how eorrelation lengths inerease with 
inereasing doping, whereas the type of anisotropy evolves from being biaxial for p < 0.12 to 
being uniaxial for p > 0.12. In order to establish a link to the domain meso-strueture and show 
that the experimentally-found peak profiles are ineompatible with mieroseopie bidireetional 
(eheekerboard) order, we define the eorrelation lengths along the x and y (equivalently H and 
K) direetions as and respeetively (the superseript refers to the peak at or Q^). The 
latter quantities are related to the longitudinal and transverse eorrelations via the relations: 

^d _ ^d 

>11 

^l=^S (Sii) 

Therefore, for domains with waveveetor = (0.31, 0) [Q;, = (0, 0.31)], .^n and are evaluated 
at a = 0° [90°] and a = 90° [180°] respeetively. For all the investigated eharge ordering peaks, 
we find that A iy Whenever > iy Hy > ix), we will talk of x-oriented (y-oriented) 
domains, and we will similarly say that the anisotropy is along x (y). Sinee we are dealing 
with local, not global eorrelations, two situations may arise: (i) a biaxial anisotropy, where 
both X- and y-oriented domains are present; (ii) a uniaxial anisotropy, where only x-oriented 
(or y-oriented) are found. Case (ii) would effeetively eorrespond to a global, maeroseopie 
anisotropy. These two possibilities are pietorially illustrated in Figs. S|7]4-D. In partieular, the 
uniaxial anisotropy shown in Figs.S|^,D assumes the presenee of y-oriented domains. Note 
that these domains need not neeessarily lie in the very same layer, but they need to be present 
at the same time within the bulk of the material (i.e., they ean be present in alternating layer, 
for instanee). The momentum-spaee representation of the order parameter - and therefore of 
the eleetronie density fluetuations is shown in the eorresponding panels in Fig. S|^-H. Here 
{pa ■ Pa(r)) (red stripes) and {pb ■ Pfe(r)) (blue stripes) represent the density-density eorrelation 
funetions for eharge modulations along a and b, respeetively. The eorresponding strueture fae- 
tors in momentum spaee 5'a(Q) and S'fe(Q) are represented as ellipsoids to refleet the presenee 
of anisotropie eorrelations. 
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By inspecting the diagrams in Figs. S|^-H, one can recognize a common trait of checker¬ 
board structures in momentum space, in that the following conditions have to be always verified 
by symmetry: 


C = C'' 

^x 

ca ^ tb 
Sy Si/? 


(S12) 


i.e., the correlation lengths along a given direction have to coincide for the two structure factors 
S'a(Q) and S'f,(Q). At the bottom-left comer of the various diagrams we also report the condi¬ 
tions that apply to the correlation lengths in the various cases presented here (see also Table S[^. 
From this theoretical analysis, we can conclude that/or both uniaxial and biaxial anisotropy it 
is in principle possible to discriminate between a pure checkerboard and a pure stripy order, 
unless in which case the patterns shown in Figs. S|^ and H become indistinguishable. 

Even though the uni- vs bi-directional character can be assessed regardless of the character of 
the anisotropy, the condition of biaxial anisotropy more intuitively discriminates between stripe 
and checkerboard charge order, as a consequence of a simple fact: for stripy nanodomains, 
the direction of the intra-domain wavevector can be locked to the direction of anisotropy - the 
two can be parallel (as in Fig. S|7]4) or perpendicular (not shown). Such configuration cannot 
be realized when the underlying order is checkerboard, since each single domain possesses 
both order parameters and therefore two wavevectors; this condition forbids by construction 
any locking to the axis of anisotropy. The case of domains which are oriented obliquely with 
respect to the crystallographic axes, but still giving rise to ellipsoidal structures in momentum 
space, can be equivalently treated by redefining the principal axes (x', y') for the domain stmc- 
ture. Ultimately, the same constraints that we derived for the correlation lengths along x and y 
would apply for x' and y' and the peak structure at S'a(Q) and S'b(Q) should still bare the same 


reciprocity conditions expressed in Eqs. S12 (but now referred to the x' and y' axes). This case 
is still not verified experimentally, and therefore rules out the presence of checkerboard order 
even in the more general case of charge order domains with different orientations. 

The comparison between the conditions laid out in Table S[^ (second column) and the val¬ 
ues for the correlation lengths determined experimentally (fourth column in Table S[^ directly 
shows that the conditions for checkerboard order are violated for all YBCO samples - indicating 
the presence of stripe order in the doping range surveyed in this study - and at the same time 
allows classifying the type of charge order (see last column in Table S[^. This classification 
indicates that Y651 and Y667 exhibit a biaxial stripy character, while Y675 can be catego¬ 
rized as having a uniaxial stripy character. This difference, as well as the fact that Y675 ef¬ 
fectively displays a global anisotropy with charge order domains preferentially elongated along 
a given direction, might relate to: (i) the crystal structure itself, with the increasing role of 
the orthorhombicity for increasing hole doping (the orthorhombic ratio a/h increases mono- 
tonically); (ii) the proximity of the chain order peak, which in Y675 (Ortho-III) is located at 
Q = 0.33r.l.u.. Both aspects might lead to a stronger influence of the orthorhombicity and its 
intrinsically anisotropic, U 4 -symmetry-broken structure onto the correlations and the domain 
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structure in the Cu 02 planes. 


Domain canting and its effect on the charge order structure in reciprocal space. So far the 

charge order has been assumed to be strictly on-axis, i.e. a modulation of the electronic density 
with crests (and valleys) parallel to the a or h axis for the and component, respectively. 
The existence of domains with slightly canted density modulations cannot be a priori ruled out, 
hence we here analyze its effect on the momentum structure of charge order and its relationship 
to the peak shapes discussed in the previous section and their link to the inner structure of charge 
order (checkerboard vs. stripes). Figure S[^ and Sj^ show the difference between on-axis and 
canted domains (with canting angle 00 = 10° in Figure S[^). 

Subsequently, we simulate the effect of a distribution of canted domains on the structure 
factor S {QxiQx)- Such distribution is Gaussian and centered about 9c = 0° (see Fig. S[^, 
angular spread is = 10°); if otherwise, the charge order peaks would be displaced away from 
the reciprocal H and K axes, at variance with our experimental observations. Figure S[^ shows 
the structure factor corresponding to a biaxial checkerboard structure with peak widths AQ“ = 
AQy = AQI, = AQy ~ 0.033 r.l.u. (see Fig. s[^ for a graphical definition), in the absence 
of canted domains. The effect of a distribution like that of Fig. S[^ is showcased in Fig. S(^, 
while Fig. sj^ shows a simple (no canting) biaxial stripe scenario with AQ“ = AQy = 0.033 
and AQI = AQy = 0.066 r.l.u. (Figs s[^-I show the same comparison for the case of unixial 
order, which also implies anisotropic correlation lengths and peak shape). The effect of the 
canting angle distribution is immediately evident in that it broadens the original charge order 
peaks in the transverse direction, as expected since canting displaces the peak along a circle 
centered at Q = 0 and with radius Q = Qco- Also, for the same reason, the peaks tend to get 
curved at their edges, an effect which becomes particularly evident when erg >5°. 

A qualitative analysis reveals a striking similarity between panels E and F, thus raising 
the possibility that, contrary to the conclusions reached in the previous section, the pattern 
shown in Figs. S[^ might in fact be compatible with checkerboard, when a distribution of canted 
domains is present. However, a quantitative analysis of the experimental peak widths shows that 
a checkerboard scenario has to be still ruled out. In the following, we consider the more general 
checkerboard case, for which these relations must hold (see also Supplementary Table S[^: 


AQl, = AQl,^AQl=AQl, 


(S13) 


The subscript i means intrinsic, as the above relations refer to the native broadening of the 
charge order peaks resulting from finite spatial correlations. As mentioned earlier, the effect of 
a distribution of canted domains will change the transverse width. The momentum broadening 
associated to the distribution of canted domains {AQ^), in the case of small angular spreads, can 
be approximated as AQ^. = 2(5co ■ (^e which, for ae = 5° = 0.087rad, gives AQ^ ~ 0.054r.l.u.. 
The analysis of Fig. S[^ demonstrates that the extra broadening adds up in quadrature to the 
intrinsic one (0.068^ ~ 0.054^ + 0.033^), by comparing the FWHM of the CO peaks along the 
transverse direction (see inset) before (0.033 r.l.u.) and after (0.068r.l.u.) the inclusion of a 
domain distribution with ag = 5°. 
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As a result, while the experimental longitudinal widths will remain unaffected: 


= AQl, = AQl 

the transverse widths will increase, and experimentally one would observe: 


(S14) 


A<3“. = = 'J + i^Qcf = \J + (AQc)" 

AQ‘,, = \/(AQ‘,y + (AQ/ = + (^Q,f = /(AQj/Tmy (S15) 

where in the second (third) equality Eqs. |S13| (Eqs. |S14| ) have been used. By taking the square of 


Eqs. S15 it directly follows that, in the case of checkerboard order in presence of a distribution 
of canted domains centered about the crystallographic axes, the following equality must hold: 

- (A(3‘J" = (A(3‘J" - (S16) 

The set of mathematical relations introduced so far has also been summarized in Eig. Sj^. 


In addition, Eq. S16 has been evaluated for the three YBCO samples, and the results are 
reported in Supplementary Table S|^ showing that such relation is always violated. This leads 
to the conclusion that, even in the presence of a distribution of domains with canted charge 
modulations, the checkerboard scenario remains incompatible with the experimental results. 


Doping- and temperature-dependent CDW peak asymmetry vs. YBCO orthorhombicity. 

In this last section, we focus on a more detailed analysis of the CDW peak asymmetry. In this 
case, the available experimental data allow extracting the longitudinal and transverse correlation 
lengths for the YBCO dopings and CDW peaks reported in Eig. 3 of the main text. Here we ex¬ 
tend such analysis to the extraction of the peak anisotropy defined as 7 = AQu/A^j^ as a func¬ 
tion of temperature, in a range where the peak widths can be extracted with reasonable precision 
(T < 100 K). The results are shown in Eig. S|^-C and suggest a non-monotonous temperature 
dependence of ya,b, which is more evident in Y651 and Y667. In particular, the peak anisotropy 
increases upon cooling down towards T >Tc, then is maximized around Tc, and eventually re¬ 
covers below Tc, consistent with the pronounced directionality of the superconductivity-induced 
suppression of CDW correlations (see again Eig. 3 in the main text). 

The magnitude of the asymmetry, which is closely related to the degree of elongation of the 
CDW peaks in momentum space (see again Eig. 1 in the main text), allows drawing important 
considerations with respect to the role of orthorhombicity in our study. The following points 
will clarify a few key aspects of our data that directly rule out the possibility that the observed 
structures are solely driven by the structural symmetry breaking associated to the chain layer: 

• In general, in presence of isotropic correlations (measured in number of unit cells) defined 
on a square lattice, the two-dimensional structure factor is also isotropic. However, when 
the underlying structure and its unit cell become orthorhombic, the peaks will develop 
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an elongation with an associated anisotropy 7 equivalent to the orthorhombic ratio a/h 
(or h/a), which is equal to a/h = 0.9902 ih/a = 1.0098) in Y651, a/h = 0.9856 {h/a = 
1.0145) in Y667, and a/h = 0.9845 {h/a = 1.0157) in Y675. The experimental values 
shown in Fig. S|^-C in almost all cases strongly deviate from the value found above for 
the orthorhombic ratio, indicating a much larger anisotropy. Furthermore, the doping 
dependence of the anisotropy, which switches from biaxial to uniaxial across p = 0 . 12 , 
is also inconsistent with a behavior completely driven by the structural orthorhombicity 
(we note again that our samples are fully detwinned). 

• If we assume that the influence acted by the chain layer and its associated potential onto 
the CUO 2 plane is more complex and, e.g., involves native anisotropic correlations in the 
CuO chain order, than one might wonder whether this alone could explain the observed 
anisotropy in the peak profiles. If this was the case, then the CDW peak along a and 
b should exhibit the very same elongation, since the chain layer breaks C 4 in a unique 
manner. In other words, the anisotropy measured at Qq and Qf, must be the same within 
such scenario, again at variance with the experimental data presented in Fig. S|^-C. 

• The deviation of the peak anisotropy parameter 7 from the isotropic case (7 = 1) is largest 
for the most underdoped sample, and is reduced for increasing hole doping (see Fig. S|^), 
in stark contrast with the evolution of the crystal orthorhombicity which instead grows as 
optimal doping is approached. 

• The temperature evolution of the peak anisotropy shows an anomaly near the supercon¬ 
ducting transition Tc. Such a behavior in principle can not be immediately deemed to 
be incompatible with the evolution of the lattice parameters a and h, whose temperature- 
dependent espansivity also shows a cusp near Tc [36]. However, the anomaly asso¬ 
ciated to the crystal orthorhombicity [36] gets more and more pronounced for increasing 
doping (with a maximum at optimal doping reflecting the complete oxygenation of the 
chains), and is therefore inconsistent with our finding of a charge order peak anisotropy 
that is more enhanced in the underdoped region, and weakens above p = 0 . 12 . 

• In addition to these aspects, a recent study by Achkar et al. [47] directly showed that the 
CDW correlations are to a large degree unaffected upon disordering the CuO chains, and 
therefore are intrinsic to the CUO 2 planes. This once again validates the scenario where 
the ordering mechanisms in the planes are largely insensitive to the phenomenology in 
the chain layer, where the orthorhombicity originates from. 

Finally, Figure S|^ shows the C 4 symmetry breaking parameter as a function of tem¬ 
perature below 100 K, for Y651 and Y667 (in Y675 no T-dependent data for are available, 
and therefore /S.C 4 cannot be evaluated). The amount of C 4 symmetry breaking ranges around 
0.2 for both doping levels, although both the error bars and the data scatter are larger for Y651, 
due to the generally lower CDW peak intensity. In any case, we can conclude that does 
not show a pronounced temperature dependence, which seems to suggest that while the overall 
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amplitude of density correlations in the Cu02 planes changes with temperature, their stripe-like 
character does not instead change significantly. 
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Supplementary Figures 



Photon Energy (eV) 

Figure SI: X-ray absorption at the Cu-Ls edge in Y675. Shown are the x-ray absorption 
profiles measured in both total fluorescenee yield (TFY) and total eleetron yield (TEY) mode, 
for ineoming vertieal (a) and horizontal (vr) polarization. The arrow marks the photon energy of 
the main transition involving planar Cu atoms (maximizing the CDW signal, see Refs. 10,30), 
where all measurements were performed. 
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Figure S2: Azimuthal geometry and angular resolution. (A) Schematics of momentum- 
resolved electronic density in YBCO systems, highlighting: (i) the chain-ordering feature (black 
dot), which in detwinned samples is only located along the reciprocal H axis; (ii) plane-ordering 
CDW peaks at ±Qa = (±0.31,0) along H (red diffuse circles); (iii) plane-ordering CDW 
peaks at ±Qfe = (0, ±0.31) along K (blue diffuse circles). The zero of the azimuthal angle (a) 
scale defines a generic direction in reciprocal space which is parallel to the H axis (top-right 
inset). Also shown are the Q-cuts for various values of a around the two CDW peaks Qa and 
Qb- (B) The top part of the graph reproduces Fig. la from the main text, showing the RXS 
scans (markers) across the CDW peak Q;, = (0, 0.31) as a function of K, with Lorentzian fits 
overlaid. The bottom part (red curve) shows a scan of the photodetector across the Bragg peak 
Qsragg = (0,1,3) as a function of detector angle 9det- The angular resolution, as obtained by 
fitting the red curve with a Gaussian lineshape, is equal to A^det = 0.2°. (C) RXS probing 
geometry highlighting how the surface-projected momentum resolution AQ|| is derived starting 
from the angular resolution A^det- 
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Figure S3: RXS data analysis. (A) Top part: RXS scans acquired on a Y667 sample around 
Qfc = (0,0.31,1.5) and at a = 90° using vertical (a) light polarization (black markers) and 
horizontal (vr) light polarization (grey markers). Fitted prohles based on Eq. S5 are shown 
as continuous lines, while the cubic baseline is shown in orange; double arrows denote the 
fuh-width-at-half-maximum, which is constrained to be equal for both scans. Bottom part: 
Background-subtracted scans, showing only the resonant contribution from the density modu¬ 
lations. (B) Same as (A), but acquired at an azimuthal position a = 0°. 
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Figure S4: Momentum-space evidence of rotational symmetry breaking. (A) Schematic rep¬ 
resentation of the comparison between the CDW structure in momentum (full ellipses) and its 
rotated version (hollow ellipses), clearly showing a violation of Eq. S7 Right panel: graphical 
definition of Ac '4 as the ratio between the light blue and grey areas as obtained from the original 
and rotated CDW structures shown in (a), taking as example the specific case of Y675. Ac 4 is 


formally defined as Ac 4 = S'(Qa) U S'(Q6) - S'(Qa) fi ^(Qb) / S'(Qa) U S'(Qfe) 


or equiv¬ 


alently as the area of the union of the two CDW peaks (one of which rotated) S'(Qa) U S{Qb) 
minus their overlap area S'(Qa) n^(Qfe), divided by their union. (B) Doping dependence of Ac 4 
(blue markers) and of the peak anisotropy (from the ratio between longitudinal and transverse 
peak width AQJAQj_) for both CDW peaks (dark red) and Q;, (red). 
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Figure S5: CDW peak linewidth vs. azimuthal angle. (A-C) Plots of the half-width-at- 
half-maximum AQ of the CDW peak Qa = (0.31, 0,1-5) as a funetion of azimuthal angle 
a, for Y651, Y667, and Y675, respeetively. (D-F) Same as (A-C), but for the CDW peak 
Qh = (0,0.31,1.5). Experimental points (colored markers) are expressed in A x 10“^. The 
uncertainties (error bars) on the data points have been estimated from the fitting analysis of RXS 
scans using standard procedures [48]. Continuous lines are fits to the /S.Q{a) data points using 
Eq. |S10 while shaded areas represent the 95% confidence bands, i.e. the range within which 
the model fit covers a 95% probability of representing the true model, in the presence of the 
reported uncertainties on AQ (indicated by error bars). The inset diagrams report the values of 
AQ along the reciprocal axes H and K, and illustrate the resulting peak shape in momentum 
space. 
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Figure S6: CDW correlation length vs. azimuthal angle. (A-C) Plots of the CDW eorrelation 
length ^ for Qa = (0.31, 0,1.5) as a function of azimuthal angle a, for Y651, Y667, and Y675, 
respectively. (D-F) Same as (A-C), but for Qb = (0,0.31,1.5). Reported in each panel are the 
values of longitudinal (,^||) and transverse (,^_l) correlations. The uncertainties (error bars) on 
the data points have been derived from the corresponding uncertainties in the values of AQ (see 

Fig. 


I using error propagation, i.e. 


OAQ 

AQ 
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Figure S7: Stripe and checkerboard charge order parameters in presence of anisotropy. 
(A,B) Stripe nanodomains in presenee of biaxial (A) and uniaxial (B) anisotropie eorrela- 
tion lengths. (C,D) Checkerboard nanodomains in presence of biaxial (C) and uniaxial (D) 
anisotropic correlation lengths. (E-H) Corresponding structure factors in reciprocal (Q) space. 
In the case of uniaxial anisotropy (F,H), a stripe and checkerboard order can be distinguished in 
Q-space provided that correlations along x or y are not equal between the two order parameter 
components, i.e. or In the case of biaxial anisotropy (E,G), stripe order pos¬ 

sesses an intra-domain wavevector which can be locked to the correlation length, whereas this 
is not possible for the checkerboard case; in such case, stripe and checkerboard charge order 
can be distinguished even more clearly. 
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Figure S8: Domains with canted modulations and their effect on the CDW structure factor. 
(A,B) On- and off-axis (10° canting angle) checkerboard domain, respectively. (C) Distribution 
of canted domains, used in subsequent simulations. (D-F) Simulation of the structure factor 
S {Qx, Qy) for a biaxial checkerboard-like charge order purely on-axis (D), or with a Gaussian 
distribution of canted domains having a = 10° (E); and (F) for biaxial on-axis stripe-like order. 
(G-I) Same as (D-F), but for uniaxial order. (J) Broadening induced by a distribution of canted 
domains with a = 5°, as seen via the linecut of S {Qx^ Qy) as indicated in the inset. (K) 
Definition of the peak widths for (red) and Qb (blue), from the region highlighted in (H). 
(L) Relations linking the experimental peak w^^th to the intrinsic broadening AQ^jy ^ 

and to the contribution due to the canted domains {AQJ) for the case of checkerboard order. 
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Figure S9: CDW peak anisotropy and C 4 symmetry breaking vs. temperature. (A-C) 

Peak anisotropy 7 ^ = (blue markers) and 7 ;, = AQ^y/AQ^ (red) for Y651, 

Y667, and Y675, respeetively. The uneertainties (error bars) on line width data points AQ 
have been estimated from the fitting analysis of RXS seans using standard proeedures [48], and 
subsequently eonverted to uneertainties for the anisotropy 7 using error propagation. Grey 
dashed lines are guides-to-the-eye. (D) € 4 ^ symmetry breaking order parameter A^^ as a 
funetion of temperature, for Y651 (red) and Y667 (light blue), respeetively. Uneertainties 
(error bars) on the values of A^^ have been estimated by taking the maximum margin af¬ 
ter having evaluated Ac^ at the extrema of the error range, i.e. for all eombinations of: 
{^Qa,|| ±^aq,,||;A(5„ _l AQ^ y ± ^ AQfe _L ± 
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Supplementary Tables 


Order type 

Conditions on correlations 

Conditions on peak width 

Stripy biaxial 

Ca > ca. cb < Cb 
<. Sy’ Sx S> Sy 

Ag“ ^ Ag“; AQl ^ AQl 

Stripy uniaxial 

Ca > ca. cb > cb 
^x Sy’ '^x 

AQ: ^ AQl; AQl ^ AQ^ 

Checkerboard biaxial 

ca ca ch ch 

^x Sa; Sy 

AQl = AQl = AQl = AQl 

Checkerboard uniaxial 

ca cb ^ ch ca 

^x Sx <I ^y ^y 

AQl = AQl^AQl = AQl 


Table S1: Constraints on correlation lengths and peak elongations for checkerboard and 
stripe orders. Reported in this table are the conditions that restrain the correlation lengths 
between the two directions and charge order peaks based on general considerations. 


Sample 

Doping 

Wavevector 

Correlations (A) 

Order type 

Y651 

Y651 

~0.10 

Qa 

Qb 

e. = 43(±7) >4 = 33(±1) 
e. = 27(±l) <4 = 41(±2) 

Stripy biaxial 

Y667 

Y667 

~0.12 

Qa 

Qb 

ex = 45(±2) >4 = 42(±2) 
e. = 35(±l) <4 = 51(±3) 

Stripy biaxial 

Y675 

Y675 

-0.13 

Qa 

Qb 

e. = 47(±4) <4 = 55(±4) 
ex = 38(±l) <4 = 54(±2) 

Stripy uniaxial 


Table S2: Charge order classification in underdoped YBCO compounds. Reported in this 
table are the correlation lengths as extracted from the RXS data. Comparison with the momen¬ 
tum structure of the order parameter allows classifying the different types of order (rightmost 
column). 
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Sample 

{AQl,f - {AQlf 


Y651 

32 (±25) X 10-® 

77 (±13) X 10-® A“^ 

Y667 

18 (±9) X 10"^ 

36 (±9) X 10"^ A“^ 

Y675 

-4 (±9) X 10-® A“^ 

23 (±9) X 10“^ A“^ 


Table S3: Evaluation of condition for checkerboard order in the presence of canted do¬ 
mains. In order for a checkerboard+canted domains scenario to be compatible with experimen¬ 
tal data, the numbers on the second and third column must be the same. When these expressions 
are evaluated and compared with the respective error bars, we find that the equality is violated 
for all three YBCO compounds. 
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